Constructing Metropolis-Hastings proposals using damped BFGS updates

The computation of Bayesian estimates of system parameters and functions of them on the basis of observed system performance data is a common problem within system identification. This is a previously studied issue where stochastic simulation approaches have been examined using the popular Metropolis--Hastings (MH) algorithm. This prior study has identified a recognised difficulty of tuning the {proposal distribution so that the MH method provides realisations with sufficient mixing to deliver efficient convergence. This paper proposes and empirically examines a method of tuning the proposal using ideas borrowed from the numerical optimisation literature around efficient computation of Hessians so that gradient and curvature information of the target posterior can be incorporated in the proposal.Comment: 16 pages, 2 figures. Accepted for publication in the Proceedings of the 18th IFAC Symposium on System Identification (SYSID

SLC26A9 Gene Is Associated With Lung Function Response to Ivacaftor in Patients With Cystic Fibrosis

Ivacaftor is a drug used to treat cystic fibrosis (CF) patients carrying specific gating CFTR mutations. Interpatient variability in the lung response has been shown to be partly explained by rs7512462 in the Solute Carrier Family 26 Member 9 (SLC26A9) gene. In an independent and larger cohort, we aimed to evaluate whether SLC26A9 variants contribute to the variability of the lung phenotype and if they influence the lung response to ivacaftor. We genotyped the French CF Gene Modifier Study cohort (n = 4,840) to investigate whether SLC26A9 variants were involved in the lung phenotype heterogeneity. Their influence in the response to ivacaftor was tested in the 30 treated patients who met the inclusion criteria: older than 6 years of age, percent-predicted forced expiratory volume measured in 1 s (FEV1pp) in the 3 months before treatment initiation ranging between 40 and 90%. Response to treatment was determined by the change in FEV1pp from baseline, averaged in 15–75 days, and the 1st-year post-treatment. We observed that SLC26A9 variants were not associated with lung function variability in untreated patients and that gain of lung function in patients treated with ivacaftor was similar to clinical trials. We confirmed that rs7512462 was associated with variability in ivacaftor-lung response, with a significant reduction in lung function improvement for patients with the C allele. Other SLC26A9 SNPs also contributed to the ivacaftor-response. Interindividual variability in lung response to ivacaftor is associated with SLC26A9 variants in French CF patients. Pharmacogenomics and personalized medicine will soon be part of CF patient care

Real-time Estimates in Early Detection of SARS

A statistical method can be used for early monitoring of the effect of disease control measures

Epidemiological evidence of higher susceptibility to vCJD in the young

BACKGROUND: The strikingly young age of new variant Creutzfeldt-Jacob disease (vCJD) cases remains unexplained. Age dependent susceptibility to infection has been put forward, but differential dietary exposure to contaminated food products in the UK population according to age and sex during the bovine spongiform encephalopathy (BSE) epidemic may provide a simpler explanation. METHODS: Using recently published estimates of dietary exposure in mathematical models of the epidemiology of the new variant Creutzfeldt Jacob disease (vCJD), we examine whether the age characteristics of vCJD cases may be reproduced. RESULTS: The susceptibility/exposure risk function has likely peaked in adolescents and was followed by a sharp decrease with age, evocative of the profile of exposure to bovine material consumption according to age. However, assuming that the risk of contamination was proportional to exposure, with no age dependent susceptibility, the model failed to reproduce the observed age characteristics of the vCJD cases: The predicted cumulated proportion of cases over 40 years was 48%, in strong disagreement with the observed 10%. Incorporating age dependent susceptibility led to a cumulated proportion of cases over 40 years old of 12%. CONCLUSIONS: This analysis provides evidence that differential dietary exposure alone fails to explain the pattern of age in vCJD cases. Decreasing age related susceptibility is required to reproduce the characteristics of the age distribution of vCJD cases

Use of Viremia to Evaluate the Baseline Case Fatality Ratio of Ebola Virus Disease and Inform Treatment Studies: A Retrospective Cohort Study.

BACKGROUND: The case fatality ratio (CFR) of Ebola virus disease (EVD) can vary over time and space for reasons that are not fully understood. This makes it difficult to define the baseline CFRs needed to evaluate treatments in the absence of randomized controls. Here, we investigate whether viremia in EVD patients may be used to evaluate baseline EVD CFRs. METHODS AND FINDINGS: We analyzed the laboratory and epidemiological records of patients with EVD confirmed by reverse transcription PCR hospitalized in the Conakry area, Guinea, between 1 March 2014 and 28 February 2015. We used viremia and other variables to model the CFR. Data for 699 EVD patients were analyzed. In the week following symptom onset, mean viremia remained stable, and the CFR increased with viremia, V, from 21% (95% CI 16%-27%) for low viremia (V < 104.4 copies/ml) to 53% (95% CI 44%-61%) for intermediate viremia (104.4 ≤ V < 105.2 copies/ml) and 81% (95% CI 75%-87%) for high viremia (V ≥ 105.2 copies/ml). Compared to adults (15-44 y old [y.o.]), the CFR was larger in young children (0-4 y.o.) (odds ratio [OR]: 2.44; 95% CI 1.02-5.86) and older adults (≥ 45 y.o.) (OR: 2.84; 95% CI 1.81-4.46) but lower in children (5-14 y.o.) (OR: 0.46; 95% CI 0.24-0.86). An order of magnitude increase in mean viremia in cases after July 2014 compared to those before coincided with a 14% increase in the CFR. Our findings come from a large hospital-based study in Conakry and may not be generalizable to settings with different case profiles, such as with individuals who never sought care. CONCLUSIONS: Viremia in EVD patients was a strong predictor of death that partly explained variations in CFR in the study population. This study provides baseline CFRs by viremia group, which allow appropriate adjustment when estimating efficacy in treatment studies. In randomized controlled trials, stratifying analysis on viremia groups could reduce sample size requirements by 25%. We hypothesize that monitoring the viremia of hospitalized patients may inform the ability of surveillance systems to detect EVD patients from the different severity strata

Influenza Epidemics in the United States, France, and Australia, 1972–19971

Influenza epidemics occur once a year during the winter in temperate areas. Little is known about the similarities between epidemics at different locations. We have analyzed pneumonia and influenza deaths from 1972 to 1997 in the United States, France, and Australia to examine the correlation over space and time between the three countries. We found a high correlation in both areas between France and the United States (correlation in impact, Spearman’s ρ = 0.76, p < 0.001, and test for synchrony in timing of epidemics, p < 0.001). We did not find a similar correlation between the United States and Australia or between France and Australia, when considering a systematic half-year lead or delay of influenza epidemics in Australia as compared with those in the United States or France. These results support a high correlation at the hemisphere level and suggest that the global interhemispheric circulation of epidemics follows an irregular pathway with recurrent changes in the leading hemisphere

S. pneumoniae transmission according to inclusion in conjugate vaccines: Bayesian analysis of a longitudinal follow-up in schools

BACKGROUND: Recent trends of pneumococcal colonization in the United States, following the introduction of conjugate vaccination, indicate that non-vaccine serotypes tend to replace vaccine serotypes. The eventual extent of this replacement is however unknown and depends on serotype-specific carriage and transmission characteristics. METHODS: Here, some of these characteristics were estimated for vaccine and non-vaccine serotypes from the follow-up of 4,488 schoolchildren in France in 2000. A Bayesian approach using Markov chain Monte Carlo data augmentation techniques was used for estimation. RESULTS: Vaccine and non-vaccine serotypes were found to have similar characteristics: the mean duration of carriage was 23 days (95% credible interval (CI): 21, 25 days) for vaccine serotypes and 22 days (95% CI: 20, 24 days) for non-vaccine serotypes; within a school of size 100, the Secondary Attack Rate was 1.1% (95% CI: 1.0%, 1.2%) for both vaccine and non-vaccine serotypes. CONCLUSION: This study supports that, in 3–6 years old children, no competitive advantage exists for vaccine serotypes compared to non-vaccine serotypes. This is an argument in favour of important serotype replacement. It would be important to validate the result for infants, who are known to be the main reservoir in maintaining transmission. Overall reduction in pathogenicity should also be taken into account in forecasting the future burden of pneumococcal colonization in vaccinated populations

Online detection and quantification of epidemics

<p>Abstract</p> <p>Background</p> <p>Time series data are increasingly available in health care, especially for the purpose of disease surveillance. The analysis of such data has long used periodic regression models to detect outbreaks and estimate epidemic burdens. However, implementation of the method may be difficult due to lack of statistical expertise. No dedicated tool is available to perform and guide analyses.</p> <p>Results</p> <p>We developed an online computer application allowing analysis of epidemiologic time series. The system is available online at <url>http://www.u707.jussieu.fr/periodic_regression/</url>. The data is assumed to consist of a periodic baseline level and irregularly occurring epidemics. The program allows estimating the periodic baseline level and associated upper forecast limit. The latter defines a threshold for epidemic detection. The burden of an epidemic is defined as the cumulated signal in excess of the baseline estimate. The user is guided through the necessary choices for analysis. We illustrate the usage of the online epidemic analysis tool with two examples: the retrospective detection and quantification of excess pneumonia and influenza (P&I) mortality, and the prospective surveillance of gastrointestinal disease (diarrhoea).</p> <p>Conclusion</p> <p>The online application allows easy detection of special events in an epidemiologic time series and quantification of excess mortality/morbidity as a change from baseline. It should be a valuable tool for field and public health practitioners.</p

Modelling the Effects of Population Structure on Childhood Disease: The Case of Varicella

Realistic, individual-based models based on detailed census data are increasingly used to study disease transmission. Whether the rich structure of such models improves predictions is debated. This is studied here for the spread of varicella, a childhood disease, in a realistic population of children where infection occurs in the household, at school, or in the community at large. A methodology is first presented for simulating households with births and aging. Transmission probabilities were fitted for schools and community, which reproduced the overall cumulative incidence of varicella over the age range of 0–11 years old